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X-RAY COMPUTED TOMOGRAPHY (CT) SYSTEM 
FOR DETECTING THIN OBJECTS 



Background of the Invention 
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The invention is directed to an improved technique 
for the detection of thin objects, for example, 
explosives along the walls of luggage . More 
specifically, the invention is directed to an improved 
technique for the detection of thin objects such as 
explosives using X-ray computed tomography (CT) . 

Detection of explosives in luggage is an extremely 
challenging problem because the amount of explosive 
required to do catastrophic damage is relatively small 
and because plastic explosives can be formed into almost 
any desired shape. Perhaps the most challenging 
configuration for detection is sheet explosive, where the 
material is deformed into a thin sheet with a very small 
physical extent in one direction. 

One conventional way of detecting explosives is 
through the use of X-ray CT. X-ray CT is a technique 
which determines the internal make-up of an object by 
passing X-rays through the object and measuring the 
attenuation of the X-rays passing through the object. In 
this technique the object is sub-divided into many 
voxels, a voxel being the basic volumetric unit for 
imaging purposes. Compared with other objects in 
luggage, explosives have a specific range of densities, 
for example, 1.2 to 1.8 gm/cc, and accordingly attenuate 
X-rays differently than non-explosives. 

In general, CT systems are designed so that the voxel 
size roughly corresponds to the smallest object of 
interest in the image. In cases where high contrast 
sensitivity is required, this approach is clearly 
justified. Indeed, voxel sizes somewhat smaller than the 
spatial dimension of interest are often used. However, 
this approach greatly increases system cost and 
complexity because it requires a large number of detector 
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elemen'ts, view angle positions, and voxels for image 
acquisition and reconstruction. X-ray source loading is 
also significantly increased due to the need to maintain 
roughly the same number of X-rays and hence the same 
signal-to-noise ratio for the smaller voxel dimensions. 

If the thin dimension of the sheet of explosive is 
smaller than the linear voxel dimension in a CT image, 
the measured density of a voxel of interest decreases due 
to the fact that the voxel is not completely filled with 
explosive. Figures 1 and 2 illustrate this problem for 
a configuration with an explosive density p of 1.5 gm/cc. 
Figure 1 shows a voxel completely filled with 

explosive, wherein the average density of the voxel is 
1.5 gm/cc. Figure 2 shows a voxel containing a 

section of sheet explosive, where the thickness of the 
sheet is 2 0% of the voxel linear dimension. The average 
density p in the voxel is reduced to 0.3 gm/cc. 

Conventional CT systems would compute a density for voxel 
Vj which is less than the density expected for an 
explosive and thus would not identify voxel as 
containing an explosive. 

The challenge is to discriminate such a sheet 
explosive from background material in the suitcase. 

Summary of the Invention 

An object of the invention, therefore, is to 
determine the presence or absence of an object in a 
three-dimensional space when the object is thin in one 
dimension. 

Another object of the invention is to provide an 
improved technique for detecting explosives. 

Another object of the invention is to provide a 
technique for detecting explosives which minimizes the 
number of voxels required to be processed while at the 
same time providing for the detection of sheet 
explosives, for example, along the walls of luggage. 
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Yet another object of the invention is to provide a 
low cost X-ray computed tomography system for explosives 
detection. 

A technique according to the invention ascertains the 
5 presence or absence of an object in a three-dimensional 
vol\xme, such as explosives in luggage* The technique 
according to the invention employs radiation scanning of 
the three-dimensional volume to determine a property (for 
example, density) of each of a plurality of voxels 

10 representing the three-dimensional volume and identifying 
voxels having similar values of the property to identify 
a contiguous group of voxels having the similar values. 
Then, the contiguous group of voxels is identified as 
containing the object if a characteristic of the 

15 contiguous group has a predetermined value. 

Other objects, features, and advantages of the 
invention will become apparent from the detailed 
description of the invention set forth below. 

Brief Description of the Drawings 

2 0 The invention will be described in greater detail 

below with r^^erence to the accompanying drawings, 
wherein: 



r/f 



Figure 1 ^illustrates a completely filled-in voxel 
Vj having an .average density p of 1.5 gm/cc; 
25 Figure 2 y illustrates a partially filled-in voxel 

having an average density p of 0.3 gm/cc; 



Figure J / is a perspective view of a suitable 
hardware arrangement according to the invention; 

Figure 4. is an end view of a gantry suitable for use 
3 0 in the invenxion; 

Figure 5^ illustrates a harmless bag Bi with small 
regions of al specific range of densities not connected 
together ; 

Figure § illustrates a suspicious bag Bj with a 
35 large contiguous region having a specific range of 
densities; and 



Figure 7^ is a flowchart for a technique to detect 
explosives in accordance with the invention. 

Detailed Description of the Preferred Embodiments 

The invention relies on the fact that explosives have 
5 a specific range of densities as compared to other items 
such as clothing in luggage and on the fact that many 
explosives must be physically in one piece, that is, 
contiguous, in order to detonate properly. According to 
the invention, an image processor performs connected 

10 component identification and labeling to identify regions 
of a specific range of densities and connects them into 
a single volume. The size of the single volume is then 
calculated to determine if the single volume is large 
enough to cause significant damage if the single volume 

15 contained explosives. 

Because the present invention is cheaper than a 
conventional CT system it can be used, for example, for 
initial screening. A secondary inspection, for example, 
by neutron or manual inspection, can be used to confirm 

2 0 that the suspect region is really an explosive. 

In this new technique according to the invention, a 
priori knowledge is combined with image processing and 
analysis to provide a system configuration with a larger 
voxel size. Consider, for example, a sheet of explosive 
25 2 mm thick with a density of 1.5 gm/cc. In a 

conventional CT system, a voxel size of 2 mm would be 
required and 500 detector elements would be required to 
cover aim field of view. Approximately 1000 view 
angles would be needed to reconstruct the 500 x 500 image 

3 0 necessary to maintain system spatial resolution. Such a 

system is extremely complex and expensive in comparison 
with the present invention, which allows use of a larger 
voxel size* 

Figures 3 and 4 illustrate hardware suitable for use 
3 5 in a preferred embodiment of the invention. Figure 3 is 
a perspective view which shows a bag B which is to be 
inspected. The bag B is moved along a conveyor 300 in 
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direction D toward a gantry 100. The gantry 100 contains 
an X-ray source which emits X-rays into the bag as the 
bag passes through the gantry 100. The gantry 100 also 
includes a set of X-ray detectors. The detectors detect 
5 X-rays which passed through the bag B. Information from 
the detectors is sent to a processor 2 00 to determine the 
attenuation of the X-rays as they pass through the bag. 
The processor 200 includes a contiguity identification 
module 210 and an object identification module 22 0, both 

10 of which will be described in further detail below. The 
attenuation information is used by processor 200 to 
compute a density for each voxel in a three-dimensional 
image of the bag. 

Reconstruction of a three-dimensional density map of 

15 the bag requires that the bag be viewed at various 
angles. Accordingly, either the gantry 100 must be 
rotated around the bag or the bag must be rotated. In 
the arrangement shown in Figures 3 and 4, the gantry 100 
is rotated because the contents of the bag would shift if 

2 0 the bag were rotated. 

Figure 4 shows an end view of gantry 100. In the 
Figure 4 implementation of the invention, gantry 100 is 
sized to accommodate a maximxxm bag size of 100 cm x 
50 cm. The gantry 100 includes an X-ray source 10 and a 

25 detector array 20. The source 10 is a 140 kVp, 1 to 5 kW 
source and has a 30® fan angle. The detector array 20 
has 1280 detection elements in an 8 x 160 array. Each 
detector element consists of a scintillator coupled to a 
photodiode and associated current integration 

30 electronics. The voxel size in this embodiment is 0.625 
X 0.625 x 1.25 cm. The number of slices for a typical 
bag (75 cm x 50 cm x 20 cm) is 60. The total inspection 
time for a typical bag is 8 seconds based on 250 views. 
General background on CT scanning, hardware, and 

3 5 signal processing may be found in "Computed Tomography 

Part I: Introduction and Industrial Applications," The 
Journal of The Minerals > Metals & Materials Societv . 
David C. Copley, Jeffrey W. Eberhard, and Gregory A. 
Mohr, Vol. 46, No. 1, January 1994, pp. 14-26; Principles 
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of Computerized Tomographic Imaging . Avinash C, Kak and 
Malcolm Slaney (IEEE Press 1988) ; and Image 
Reconstruction From Projections , Gabor T. Herman 
(Academic Press 1980) . The entire contents of these 
5 publications are incorporated herein by reference. 

In this new technique according to the invention, a 
large voxel size on the order of 1 cm can be employed. 
A large sheet of explosive passing through a given voxel 
would fill 20% of the voxel volume, resulting in an 

10 average density of the voxel of 0.3 gm/cc. This density 
is still sufficiently large to be discriminated against 
a background bag density of 0.2 gm/cc. Thus, voxels 
within a chosen range of densities, for example, from 
0.25 to 1.8 gm/cc are identified as potentially 

15 containing explosives. These regions which have 
densities within a specific range are illustrated in 
Figures 5 and 6 as regions Rj, Rj, R3, and R4. Figures 5 
and 6 will be discussed in more detail below. 

After the above-background regions are identified, 

2 0 a three-dimensional grey-scale connected component 

identification and labeling process joins all physically 
adjacent voxels in the selected density range. Component 
identification and labeling techniques are set forth in 
Chapter 4 of Robot Vision by Berthold Klaus Paul Horn 
25 (MIT Press 1986) . This publication is incorporated 
herein by reference. 

The procedures described in this text are modified 
for use in the present invention in that the present 
invention employs a three-dimensional grey scale 

3 0 procedure instead of a two-dimensional binary procedure. 

Instead of checking to determine if a voxel is a 1, a 
check is made to determine if its difference a with 
respect to a selected reference value is less than a 
preset threshold. If it is, the voxel is treated just 
3 5 like a 1 in the binary case. In three-dimensions, 
volumes are used instead of areas. For any given voxel, 
the neighbors are divided into a plane above, a same 
plane, and a plane below the voxel of interest. A three- 
dimensional raster scan is performed, proceeding from top 
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to bottom, while a scan within a plane is a two- 
dimensional raster scan. As in the two-dimensional case, 
a subset of the neighboring elements is used in the 
labeling scheme. All voxels in the plane above are 
5 evaluated. If any voxel in the plane above is within the 
threshold difference value a of the target voxel, it is 
assigned the label of that voxel. In the plane of the 
voxel of interest, the same voxels as in the two- 
dimensional case are utilized. As in the two-dimensional 

10 case, it is possible that two different labels have been 
used for parts of one component. Indeed, this is the 
case if the two voxels which give labels to the central 
voxel of interest are connected only point-wise to the 
central voxel (not on a face or edge) . In such a case, 

15 it is necessary to note that the two labels are 
equivalent and use either of them for the central voxel. 
Re-labeling using a second scan over the image may be 
necessary. 

Voxels in the specified density range which are 

2 0 physically adjacent in three-dimensions and have a 

density variation less than a predetermined threshold are 
grouped and assigned with a label for identification. 
Because this adjacency check is performed in three- 
dimensions, thin regions of any shape in any orientation 
25 are easily identified. 

Next, the number of voxels in each region is 
determined and compared to a threshold. Small regions, 
that is, regions containing only a small number of voxels 
are rejected as being "harmless." This situation is 

3 0 illustrated in Figure 5, which shows a "harmless" bag Bi 

with small regions Ri, Rj, and Rj in the specified density 
range not connected together. Large contiguous regions, 
that is, regions containing more voxels than a preset 
threshold, are identified as suspect. This situation is 
35 illustrated in Figure 6, which shows a suspicious bag B2 
with a large contiguous region R4 of material in the 
specified density range. The mass contained in any 
suspect region (s) is then calculated by multiplying the 
volume of each voxel in the region by its density. If 




• 



-8- 



10 



15 



20 



25 



30 



35 




the resulting mass is greater than a preset threshold, 
for example, 1000 gm, the region is tentatively 



performed, either by a second inspection technique, such 
as pulsed fast neutron analysis or by opening the bag. 

A detailed example of the above-described technique 
will now be described with reference to Figure 7. Most 
of the steps shown in Figure 7 are performed in the 
contiguity identification module 210 and the object 
identification module 22 0 of processor 200. Depending on 
the specific application at hand, these modules can be 
implemented by software, hardware, or a combination of 
both. 

It should be noted that the technique illustrated in 
Figure 7 is only one example of an application of the 
principles of the invention. Those skilled in the field 
will appreciate that numerous modifications and 
variations of the Figure 7 technique are possible. 

The Figure 7 example is based on the following a 
priori information : 

(1) The explosive density is in the range of 1.2 to 
1.8 gm/cc; 

(2) The background density in the suitcase is 
approximately 0 . 2 gm/cc ; 

(3) The detection of small amounts of explosives, 
for example, 100 cc or 150 gm, is not desired; 
and 

(4) Partial volvime artifacts reduce density 
contrast in direct proportion to the lack of 
filling of a voxel. 

Suitable modifications can be made to this a priori 
information based on the specifics of the detection 
problem at hand. In the Figure 7 example, 0.3 to 
1.8 gm/cc is selected as the density range of interest. 

In step Sj the suitcase is scanned to determine a 
linear attenuation coefficient for each voxel, which in 
turn represents the density of each voxel. Suitable 
techniques for scanning and density determination are 
described in the above-cited references concerning CT 



identified as explosive. 
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scanning. This density information is stored in a memory 
for processing in steps S2 to S-j, which together identify 
potentially explosive regions in the CT image of the 
suitcase, 

5 In step S2, voxels having a similar range of 

densities, that is, similar densities between 0.3 and 
1.8 gm/cc are connected and labelled as a contiguous 
region. In step S3 the number of voxels in each 
contiguous region are counted. In step S4 the volume of 

10 each contiguous, or connected and labeled, region is 
determined by multiplying the number of voxels for that 
region (from step S3) by the voxel volume. Step S5 
determines, for each contiguous region, whether the 
volume of the contiguous region is greater than a 

15 threshold Tj, for example, 750 cc. Any region whose 
volume is greater than threshold Ti is considered suspect 
and processing continues to step Sg. 

In step S5 the mass of each suspect contiguous region 
is determined by summing over the suspect contiguous 

2 0 region the product of each voxel density and voxel 

volume. Step S7 determines whether the mass of each 
suspect contiguous region (from step S5) is greater than 
threshold T2, for example, 1000 gm. If the mass of a 
suspect contiguous region is greater than threshold T2, 
25 then the region is tentatively identified as explosive 
and processing proceeds to step Sg. Step S7 could include 
activating an alarm. In step Sg the presence, or 
absence, of an explosive is verified by an additional 
inspection method, for example, by pulsed fast neutron 

3 0 analysis or by opening the bag. 

Because this new technicjue allows the use of large 
voxels instead of 2 mm voxels, the number of detector 
elements required is reduced. For example, if 1 cm 
voxels are used, the number of detector elements required 

35 is reduced to 100, the number of view angles is reduced 
to approximately 200, and the image size is reduced to 
100 X 100. The input data set size, which is 

proportional to the niimber of detectors times the number 

^ of view angles, is therefore reduced by a factor of 25 
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(or more, because slice 'thickness can also be increased) . 
The image reconstruction time, which is proportional to 
the nxmber of view angles times the number of voxels in 
the image, is reduced by a factor of 125. These drastic 
5 reductions in data and computational load make practical 
and reliable CT systems for baggage inspection much 
simpler. 

Although the invention has been described above with 
respect to certain specific applications and 

10 implementations of the invention, the scope of the 
invention is not limited to the specific applications and 
implementations described above. Various modifications, 
variations and applications within the spirit and scope 
of the invention will occur to those skilled in the field 

15 after receiving the above teachings. For example, the 
invention is not limited to the physical arrangement 
illustrated in Figures 3 and 4. Although the invention 
is particularly useful to detect sheet explosives in 
luggage, it is generally useful whenever an object of 

2 0 interest is smaller in one direction than the linear 
dimension of the voxels in the CT image, and can be used, 
for example, to detect and characterize delaminations in 
composite materials. Accordingly, the scope of the 
invention is defined by the following claims. 
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